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Sphere  impact  experiments  are  used  to  calibrate  and  validate  ceramic  models  that  include  statistical  variability  and/ or  scale 
effects  in  strength  and  toughness  parameters.  These  dynamic  experiments  supplement  traditional  characterization  experiments 
such  as  tension,  triaxial  compression,  Brazilian,  and  plate  impact,  which  are  commonly  used  for  ceramic  model  calibration. 
The  fractured  ceramic  specimens  are  analyzed  using  sectioning,  X-ray  computed  tomography,  microscopy,  and  other  tech¬ 
niques.  These  experimental  observations  indicate  that  a  predictive  material  model  must  incorporate  a  standard  deviation  in 
strength  that  varies  with  the  nature  of  the  loading.  Methods  of  using  the  spherical  indentation  data  to  calibrate  a  statistical 
damage  model  are  presented  in  which  it  is  assumed  that  variability  in  strength  is  tied  to  microscale  stress  concentrations 
associated  with  microscale  heterogeneity. 
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Numerical  modeling  tools  have  long  been  used 
to  capture  the  behavior  of  ceramics  in  armor  systems. 
For  example,  seminal  constitutive  modeling  work  of 
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Johnson  and  Holmquist1  has  established  an  essential 
foundation  for  describing  ballistic  failure  of  ceramics. 
The  present  work  aims  to  extend  this  class  of  model  by 
incorporating  the  effect  of  microscale  heterogeneity  that 
leads  to  pressure  dependence  of  failure  mechanisms. 
Without  such  enhancements,  a  ceramic  model  that  is 
tuned  to  match  one  experiment  often  fails  to  predict  the 
results  in  a  different  application.2 

The  ceramic  experiments  of  Hauver  et  al.?  Lundberg 
and  Lundberg,4  Orphal  and  Franzen,5  and  Anderson 
et  al?  provide  penetration  rates  and  other  time-varying 
target  failure  information.  These  time-resolved  experi¬ 
ments  are  used  to  calibrate  the  damage  model  parameters 
(such  as  strain-to-failure  or  failure  energy)  that  are  asso¬ 
ciated  with  the  progression  of  failure  from  its  initiation  to 
complete  loss  in  strength.7  Failure  initiation  itself  begins 
at  statistically  distributed  weak  points  in  the  material,  thus 
requiring  statistical  spatial  variability  in  strength  that  is 
proposed  here  to  be  parameterized  via  spherical  indenta¬ 
tion  experiments.  Specifically,  the  cracking  patterns  ob¬ 
served  in  spherical  indentation  data  indirectly  quantify 
microheterogeneity. 

The  evolution  of  damage  in  ceramics  due  to 
projectile  impact  has  been  studied  extensively  in  the 
past.  Excerpts  from  ceramic  society  proceedings  such  as 
Sagamore8  and  PacRim,9  contain  extensive  references. 
Normandia  and  Gooch10  reviewed  earlier  attempts  to 
screen  ceramics  using  sphere  impact  tests  by  Donaldson 
et  al .n  and  Donaldson  and  McDonough12  and 
critiqued  by  Sternberg.13  Low-velocity  impacts  using 
tungsten  carbide  (WC)  spheres  to  generate  damage  were 
presented  by  Shockey  et  al.lA  and  Kim  et  al 15  for 
ceramic  targets.  Historical  references  of  high-velocity 
penetration  using  WC  data  appear  in  Herrmann  and 
Jones,16  Goldsmith,17  and  Martineau  et  al ,18  There  are 
a  number  of  helpful  reviews  by  Lawn,19  Stronge,20  and 
Fischer-Cripps  et  al.21  For  damage  analyses,  there  is  a 
useful  overview  of  ceramic  impacts  by  Tanabe  et  al?2 
The  current  paper  extends  the  previous  ceramic  failure 
analysis  with  modeling  techniques  that  incorporate 
experimentally  observed  variability  and  size  effects.23 

The  indentation  experiments  presented  in  this 
paper  examine  the  ceramic  damage  morphology  of  sim¬ 
plified  low-velocity  projectile  impacts.  WC  spheres  with 
velocities  between  50  and  500  m/s  impact  ceramic  tar¬ 
gets  in  metal  cups.  Ceramic  targets  consist  of  silicon 
carbide  (SiC)  and  boron  carbide  (B4C).  These  generic 
screening  experiments  allow  for  a  detailed  analysis  of 
specific  materials  of  interest,  with  the  overall  goal  being 


a  connection  of  the  ballistic  performance  to  statistically 
varying  material  properties. 

The  sphere  impact  tests  are  used  to  calibrate  a  new 
ceramic  model,  called  Kayenta,  currently  implemented 
in  the  shock-physics  finite-element  code  ALEGRA.24 
This  model  is  an  extension  of  the  Sandia  GeoModel25 
that  was  originally  developed  for  geological  materials, 
but  has  been  enhanced  to  include  scale-dependent  sta¬ 
tistical  damage  for  ceramics  and  metals.26  Scale  depen¬ 
dence  and  statistical  variability  are  included  not  only  to 
reduce  mesh  sensitivity  but  also  to  improve  predictive 
capability  by  providing  a  distribution  of  weak  points 
that  give  rise  to  localized  cracking.27 

The  Kayenta  model  in  ALEGRA  uses  a  new  failure 
distribution  methodology  that  shares  similarities  to  clas¬ 
sical  Weibull28  theory  in  hydrostatic  tension  but  which 
smoothly  deviates  from  Weibull  theory  for  other  load¬ 
ing  paths  that  involve  shear.  A  Weibull  distribution  is 
used  merely  as  an  option  that  can  be  replaced  with  other 
distributions  if  necessary.  This  approach  uses  not  only 
the  material  properties  measured  from  laboratory  ex¬ 
periments  but  also  the  length  scales  associated  with  the 
particular  experiments.  Volume  scaling  concepts  similar 
to  the  work  of  Bazant  and  Planas29  Bazant30  and  Trust- 
urn  and  Jayatilaka31  are  implemented  at  the  finite  ele¬ 
ment  level. 

Accounting  for  the  size  effects  of  strength  helps 
mitigate  mesh  dependency  problems.  Moreover,  incor¬ 
porating  spatially  variable  strengths  automatically  results 
in  an  “ event  horizon”  effect  in  which  increasing  the 
strain  rate  causes  an  increase  in  the  number  of  fractures 
because  of  the  finite  time  required  for  release  waves  to 
propagate  through  the  material.32 

To  parameterize  the  new  material  model,  the  dy¬ 
namic  characterization  work  of  Grady,33  Vogler  et  al.?A 
and  Dandekar35  is  used.  Additionally,  the  compilations 
of  Holmquist  et  al?G  are  an  essential  starting  point.  The 
original  high  pressure  yield  surface  calibration  is  based 
on  the  Johnson-Holmquist  (JH)  ceramic  models.37  The 
Kayenta  model’s  completely  failed  material  behavior  in¬ 
volves  concepts  similar  to  the  work  of  Duncan38  and 
Anderson  and  Chocron.39  Traditional  quasi-static  tests 
conducted  by  A.  A.  Wereszczak  (personal  communica¬ 
tion)  and  others  provide  additional  invaluable  informa¬ 
tion.  The  failure  distributions  used  in  the  Kayenta 
model  stem  from  the  actual  measured  uncertainty  in 
the  characterization  experiments.  New  experiments  were 
completed  at  the  Sandia  National  Labs40  to  obtain  mea¬ 
surements  of  the  failure  surface  under  triaxial  loading 
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conditions,  an  area  of  the  failure  surface  where  no  data 
were  available. 

Methods 

The  Kayenta  Model 

Kayenta  strikes  a  balance  between  first-principals  mi¬ 
cromechanics  and  phenomenological,  homogenized,  and 
semiempirical  modeling  strategies.41  A  set  of  standard 
strength  experiments  (e.g.,  plate  impact,  triaxial  compres¬ 
sion  [TXC],  diametral  compression)  is  used  to  set  model 
parameters.  Uncertainty  observed  in  these  experiments,  as 
well  as  variability  in  crack  patterns  for  spherical  indenta¬ 
tion,  are  used  to  set  values  for  the  scale  effects  and  statistical 
parameters  (a  complete  derivation  and  analysis  of  this  new 
method  is  presented  in  Brannon  and  Strack42).  Once  these 
parameters  have  been  determined,  ballistic  experiments  are 
used  to  calibrate  additional  parameters  in  the  model. 

Like  the  JH  ceramic  model,  as  well  as  other  classical 
damage  models,43  Kayenta  allows  the  stress  to  reach  a 
pressure-dependent  limit  threshold  surface  that  then 
collapses  with  damage  down  to  a  completely  failed 
strength  surface  that  is  similar  to  that  of  sand  or  gravel. 
Unlike  the  JH  model,  which  imposes  a  fixed  strain- 
to-failure,  Kayenta  imposes  a  scale-dependent  time  to 
failure  that  indirectly  regulates  the  energy  of  failure. 
Kayenta  also  includes  the  third-invariant  dependence  as 
well  as  porosity  effects  in  its  yield  surface;  see  Fig.  1.  The 
Kayenta  limit  surface  in  TXC  specifies  shear  strength  x 
dependence  on  pressure  p  according  to 

x  =  Y[\-e-s{t+h)IY]  (1) 


where  the  material  parameters  are  the  high-pressure 
limit  strength  Y,  the  hydrostatic  tensile  strength  h,  and 
the  initial  slope  5  at  the  hydrostatic  tensile  limit.  Any  of 
these  parameters  can  be  varied  statistically  and  still  gen¬ 
erate  a  convex  surface.  The  addition  of  uncertainty  to 
the  yield  function  is  another  distinguishing  feature  of 
the  Kayenta  model.  Kayenta’s  spatially  variable  proba¬ 
bility  functions  incorporate  uncertainty  in  flaw  size  and 
orientation  of  the  macroscopically  homogeneous  mate¬ 
rial.  Through  proper  scaling,  the  behavior  is  captured  in 
the  fuzzy  failure  boundary  theory  that  incorporates 
length  scales  to  mitigate  mesh  dependence.44  The  fail¬ 
ure  distributions  predict  that  larger  volumes  are  inher¬ 
ently  weaker  and  fail  earlier  than  smaller  ones  due  to  the 
probability  of  critical  flaws  (Fig.  2). 

The  Weibull  distribution  is  applied  to  the  horizon¬ 
tal  pressure  axis  of  Fig.  2  with  h,  which  is  the  peak 
strength  in  hydrostatic  tension.  The  parameter  to  be 
perturbed,  in  this  case  the  initial  peak  h  value,  is  im¬ 
plemented  at  the  finite  element  level  by  the  following 
equation: 


Vln  R  11A 

Kln(l/2) 


where  h  is  the  median  value  of  the  property  h  as  mea¬ 
sured  using  a  sample  volume  V,  R  is  a  random  number, 
and  m  is  the  Weibull  modulus  quantifying  the  observed 
variability  in  the  material  property.  The  Weibull  distri¬ 
bution  is  seeded  about  the  median  value  for  the  mate¬ 
rial,  with  the  value  adjusted  for  the  finite  element 
volume  V.  A  random  distribution  of  R  values  is  used 
for  the  elements,  with  the  overall  average  being  a  factor 


Fig.  1.  (a)  Yield  surface  with  a  high-pressure  cap  shown  as  a  wireframe,  (b)  Meridional  side  view,  with  the  hydrostat  as  abscissa,  (c) 

Noncircular  octahedral  profile  allowing  for  well-known  extension— compression  differences  in  brittle  strength. 
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Fig.  2.  The  Weihull  scale  effect  is  applied  only  on  the  pressure  intercept.  The  curved  lines  on  the  left  plot  are  median  limit  surfaces  for  different 
specimen  volumes  showing  that,  as  shown  in  the  plot  on  the  right,  strength  remains  bounded for  any  stress  path,  such  as  uniaxial  stress,  having  a 
deviatoric  component. 


of  the  median,  the  finite  element  volumes,  and  the  Wei- 
bull  modulus.  As  seen  in  Fig.  2,  while  the  hydrostatic 
tension  h  can  increase  without  bounds,  the  values  of  T 
approaches  the  limit  value  of  Y.  In  this  way,  the  model 
implementation  incorporates  observed  behavior  seen  in 
the  characterization  experiments. 

Experimental  Parameterization 

A  traditional  yield  model  presumes  that  shear  strength 
is  a  deterministic  scale-insensitive  function  of  pressure.  The 
Kayenta  model  incorporates  observed  variability  in  strength 
so  that  the  failure  threshold  is  replaced  by  an  uncertainty 
band  centered  about  a  scale-sensitive  median  curve. 

As  illustrated  in  Fig.  3,  the  observed  degree  of  vari¬ 
ability  decreases  with  increasing  pressure,  thus  giving  a 


(GPa) 


Fig.  3.  Meridional  view  of  the  Kayenta  material  model  and 
experiments  for  silicon  carbide. 


thick  uncertainty  band  at  low  pressure  that  becomes  in¬ 
creasingly  thin  with  increased  confinement.  For  exam¬ 
ple,  in  the  high-pressure  plate  impact  experiments,  there 
is  very  little  variability,  corresponding  to  a  large  Weibull 
modulus.34  In  low-pressure  Brazilian  tests  performed  at 
Sandia  (see  Fig.  4),  the  Weibull  modulus  is  8.3.  Figure  3 
also  incorporates  normalization  of  the  data  to  a  com¬ 
mon  reference  scale.  Namely,  if  strength  measurements 
are  performed  using  a  specimen  volume  V,  then  refer¬ 
ring  to  Eq.  (2),  Weibull  theory  predicts  that  strength 
measured  using  a  reference  sample  volume  V  is  different 
by  a  factor  of  (V /V)l^m .  Raw  experimental  values  for  a 
variety  of  specimen  sizes  are  displayed  as  small  blue  dots 
in  Fig.  3,  and  larger  red  dots  show  the  same  data  scaled 
to  a  common  reference. 

To  date,  the  primary  experiments  validating  a  scale 
effect  in  SiC  have  been  diametral  compression,  or  Bra¬ 
zilian  experiments.  Cylinders  ranging  between  3.175 


Fig.  4.  Brazilian  Weibull  distribution  and  volume  scaling  of 
silicon  carbide. 
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and  25.4  mm  in  diameter  are  quasi-statically  com¬ 
pressed  between  two  platens  until  failure. 

Figure  4  illustrates  the  importance  of  uncertainty 
and  scale  effects.  The  abscissa  is  Brazilian  failure 
strength  on  a  logarithmic  scale.  The  ordinate  is  a  dou¬ 
ble-logarithmic  transformation  of  the  probability  that 
the  sample  is  safe  from  failure  The  solid  filled 

measured  data  points  are  fitted  to  solid  straight  lines  in 
which  the  slope  equals  the  Weibull  modulus.  One  col¬ 
ored  line  and  symbol  set  was  drawn  for  each  of  the  three 
different  sample  diameters  tested.  The  corresponding 
hollow  data  points  and  dashed  lines  are  simulations  us¬ 
ing  the  Kayenta  damage  model  with  statistics  and  scale 
effects.  The  simulations  predict  the  right  trend  of  an 
increase  in  strength  and  a  decrease  in  slope  as  the  spec¬ 
imen  size  decreases.  Upon  inspection  of  the  simulation 
results,  it  is  clear  that  they  are  highly  mesh  sensitive. 
Efforts  to  resolve  this  numerical  verification  issue  have 
included  alternative  volume  scaling  formulations,  and 
work  is  underway  to  investigate  its  ties  to  energy  release 
rates  as  well  as  under-resolution  of  stress  associated  with 
low  order  shape  functions  in  the  finite  element  model. 

Triaxial  compression  and  extension  tests  provide 
additional  information  about  the  yield  surface  and  vari¬ 
ability  of  brittle  materials.  Somewhat  a  misnomer,  “tri- 
axial”  refers  to  an  axisymmetric  stress  state  having  two 
equal  eigenvalues.  For  TXC,  the  axial  stress  (single 
eigenvalue)  is  more  compressive  than  the  lateral  stresses 
(dual  eigenvalues).  For  triaxial  extension  (TXE),  the  ax¬ 
ial  stress  is  less  compressive  than  the  lateral  stresses.  In 
most  cases,  all  the  stresses  in  TXE  are  compressive,  de¬ 
spite  the  use  of  the  term  c ‘extension”  in  the  name  of  the 
test.  Strength  in  unconfined  uniaxial  tension  is  orders  of 
magnitude  smaller  than  in  TXC.  This  effect  can  be 
modeled  purely  through  the  pressure  dependence  of 
strength.  Fiowever,  at  low  pressures,  the  strength  of  rock 
is  strongly  dependent  on  the  third  stress  invariant.  Spe¬ 
cifically,  when  comparing  TXE  and  TXC  at  the  same 
pressure ,  data  show  that  the  equivalent  shear  stress  at 
failure  in  TXC  can  be  as  much  as  twice  as  large  as  that  in 


TXE.  This  observation  requires  limit  surface  shapes  that 
have  strong  triangular  octahedral  profiles  at  low  pressure 
as  illustrated  in  Fig.  1.  Consistent  with  observations  for 
quasi-brittle  materials,  the  Kayenta  model  allows  the 
octahedral  profile  to  smoothly  transition  from  a  triangle 
(corresponding  to  a  maximum  principal  stress  or  strain 
criterion)  at  low  pressure  to  a  circle  (corresponding 
to  a  von  Mises  or  Drucker-Prager  yield  model)  at 
high  pressure. 

Plate  impact  experiments  probe  the  high-pressure 
regime  of  the  model.  As  line  Visar  experiments,  in 
which  multiple  points  of  data  are  recorded,  become 
more  prevalent,45  more  statistical  information  can  be 
discerned  from  a  single  test.  The  current  series  of  ex¬ 
periments  by  Vogler  et  al ,34  supplement  strength  and 
pressure  data  consistent  with  previous  work.  These  ex¬ 
periments  provide  controlled  model  parameterization 
data  in  a  regime  that  exists  in  many  dynamic  ballistic 
experiments. 

Ballistic  Experiments 

In  this  study,  ballistic  impact  experiments  up  to 
500  m/s  are  conducted  on  SiC  and  B4C  ceramics  using 
WC  spheres.  The  test  procedure  consists  of  firing  6.35- 
mm  diameter  spheres  launched  in  a  sabot  from  a  gas 
gun.  The  WC  spheres  contain  6%  cobalt  by  volume, 
with  a  nominal  density  of  14.90  g/cm3.  Time  is  mea¬ 
sured  between  two  lasers  at  the  end  of  the  barrel  to  de¬ 
termine  velocity,  or  with  two  orthogonal  pairs  of 
150keV  X-ray  radiograph  systems.  The  target  consists 
of  a  25.4-mm  diameter  and  length  ceramic  cylinder, 
manufactured  by  BAE  Systems  (Vista,  CA).  The  two 
materials  tested  are  SiC  (SiC-N)  and  B4C.  Table  I  has  a 
summary  of  the  material  properties  for  the  target  materials. 

The  ceramic  is  slip-fit  inside  a  6.35-mm  walled 
titanium  alloy  cup,  which  is  held  in  place  by  a  rigid 
structure.  The  opening  of  the  cup  allows  the  sphere  imp- 
actor  to  strike  the  ceramic  face  directly.  The  cup  keeps 
the  damaged  ceramic  in  place  for  posttest  examination. 


Table  I.  Ceramic  Target  Properties 


Grain 

Young’s 

Knoop 

Fracture 

Flexure 

Density 

Size 

Modulus 

Poisson’s 

Hardness 

Toughness 

Strength 

Material 

(g/cm3) 

(gm) 

(GPa) 

Ratio 

(GPa) 

(MPam1/2) 

(MPa) 

Polytype 

SiC-N 

3.22 

3.3 

452 

0.165 

18.6 

5.1 

620 

6H,  15R 

b4c 

2.52 

15.0 

465 

0.170 

18.9 

2.9 

450 

B4C,  c 
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A  significant  benefit  of  this  technique  is  material  recov¬ 
ery  for  microscopic,  X-ray  computed  tomography,  and 
cross  section  analyses.  These  experiments  can  be  used  to 
calibrate  the  Kayenta  model’s  damage  evolution  and 
variability  parameters. 

Figure  5  illustrates  the  relationship  between  the 
model  implementation  of  variability,  and  an  actual  sim¬ 
ulation  of  the  ballistic  impact  experiments.  Lagrangian 
finite  elements  in  the  ceramic  target  on  the  right  are 
colored  by  the  seeded  tensile  strength  h.  The  red  color 
corresponds  to  Weibull  seeded  elements  with  the  high¬ 
est  tensile  strength  (the  largest  negative  values  on  the 
pressure  axis  in  Fig.  5).  These  elements  have  the  greatest 
strength  based  on  the  Weibull  deviation  from  the 
median  (green)  values.  As  can  be  seen  in  the  figure, 
the  target  appears  to  be  mostly  green,  corresponding 
to  the  median  strength  of  the  material.  Blue  elements 
in  the  target  correspond  to  weaker  elements,  which  used 
the  lower  blue  curve  of  the  model. 

In  addition,  the  radial  mesh  used  in  this  example  has 
slightly  smaller  elements  in  the  center  of  the  simulation, 
with  larger  elements  near  the  outer  edge  of  the  target.  The 
volume  scaling  of  the  material  strength  in  these  elements 
is  witnessed  by  the  greater  accumulation  of  higher 
strength  (yellow  and  red)  in  the  middle  of  the  simula¬ 
tion.  By  relating  the  element  size  and  all  model  calibra¬ 
tion  to  a  common  reference  volume,  the  simulation  will 
accurately  predict  the  material  performance  differences  in 
a  simulation  that  models  a  sample  with  1  mm3  elements 
versus  one  with  cubic  meter  wide  elements. 

Incorporation  of  statistics  and  scale  effects  signifi¬ 
cantly  mitigates  mesh  size  sensitivity  in  dynamic  indenta¬ 
tion  problems,  which  have  been  reported  to  have  severe 
mesh  texture  bias  in  glass  indentation  simulations.46 

Figure  6  shows  the  same  500  m/s  impact  simulation 
of  SiC  run  with  three  different  mesh  sizes.  Damage 


Fig.  6.  Mitigating  mesh  dependency  in  three  dijferent  resolution 
silicon  carbide  indentation  simulations.  Larger  to  smaller  mesh  sizes 
shown  from  left  to  right. 


(light  colored  elements)  occurs  in  roughly  the  same 
locations  and  at  the  same  time  in  all  three  mesh  reso¬ 
lutions.  However,  a  similar  reduction  in  mesh  sensitivity 
was  not  observed  in  Brazilian  testing,  even  at  high  rates. 
The  contrast  between  this  mixture  of  positive  and  neg¬ 
ative  results  indicates  that  proposed  methods  for  limit¬ 
ing  mesh  sensitivity  should  be  demonstrated  to  be 
effective  under  a  variety  of  loading  scenarios  using  a 
single  material  parameter  set. 


Results  and  Discussion 

The  ballistic  impact  experimental  calibration  and 
validation  involves  recording  the  failure  morphology  of 
a  variety  of  simplified  ballistic  impact  events.  Figure  7  is 
a  labeled  overlay  of  the  experimental  results  from  a  cross 
sectioned  SiC  target  impacted  at  375  m/s. 


Fig.  7.  Crosssectioned  silicon  carbide  specimen  with  highlighted 
features. 
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Fig.  8.  Silicon  carbide  cross  section  of 500  m/s  impact  (left)  and  zoom  in  of  comminuted  zone  (right). 


Lateral  cracks  in  Fig.  7  branch  out  horizontally 
from  the  cone  cracks  that  originate  at  the  impact  point. 
A  ring  crack  is  essentially  a  circular  crack  that,  along 
with  radial  cracks,  is  visible  from  the  top  surface.  In 
addition,  for  SiC,  a  visual  damage  region  forms  between 
the  major  cone  cracks  directly  under  the  impact  point. 
This  comminuted  region  consists  of  a  dense  population 
of  microcracking  between  the  grain  boundaries,  as  well 
as  grain  pullout  as  seen  in  Fig.  8. 

As  observed  in  the  experiments,  the  comminuted 
region  is  less  noticeable  in  the  extremely  low  velocity 
impacts.  Moreover,  the  existence  of  a  comminuted  re¬ 
gion  appears  to  be  material  dependent.47,48 

The  cone  cracking  distributions  and  angles  are  also 
of  interest  in  determining  how  the  material  will  fail  dy¬ 
namically.  As  the  impact  velocity  increases,  the  angle 
(measured  from  the  horizontal)  and  number  of  cone 
cracks  increases  substantially.  If  the  thickness  of  the  ce¬ 
ramic  is  decreased  or  if  the  projectile  impulse  is  in¬ 
creased,  the  cone  cracks  eventually  become  vertical. 


Several  samples  have  also  been  analyzed  using  X-ray 
computed  tomography.49  The  potential  to  analyze  the 
results  in  a  nondestructive  manner  and  to  generate  3D 
image  maps  of  experiments  makes  this  method  very  at¬ 
tractive.  It  is  possible,  for  example,  to  measure  the  angle 
of  the  major  cone  crack  in  one  385  m/s  impact  for  all 
360°,  which  is  approximately  66°  from  vertical. 

Macroscopically,  the  radial  cracking  of  specimens  is 
another  failure  mode  of  interest  in  determining  how  a 
ceramic  will  perform  in  a  ballistic  event.  Radial  cracks 
serve  as  symmetry  breakers,  which  can  be  realistically 
modeled  by  the  introduction  of  variability,  volume  scal¬ 
ing,  and  by  including  the  low-pressure  strength  differ¬ 
ences  discussed  earlier.  Figure  9  is  an  example  of  the 
radial  cracking  seen  in  experiments  and  simulations  with 
SiC.  Without  the  incorporation  of  variability  in  the 
simulations,  the  ceramic  material  model  would  behave 
the  same.  As  such,  when  failure  occurs,  it  catastrophi¬ 
cally  fails  everywhere  at  once,  effectively  removing  the 
brittle  failure  capabilities  that  may  be  incorporated  in  a 


Fig.  9.  X-ray  computed  tomography  image  of  top  surface  cracking  (left)  and  corresponding  simulation  ( right)  from  a  385  m/s  impact  of  silicon 
carbide. 
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Fig.  10.  Increase  in  radial  cracking  in  boron  carbide  experiments  due  to  an  increase  in  impact  velocity  from  100  m/s  (left)  to  400  mis  ( right). 


material  model.  By  including  variability,  the  numerical 
issues  that  hinder  most  simulations  are  mitigated  and 
the  model  can  react  as  designed. 

As  seen  in  the  experiments  in  Fig.  10,  the  number 
of  radial  cracks  increases  with  an  increase  in  impact  ve¬ 
locity.  The  Kayenta  model  accurately  predicts  an  in¬ 
creased  frequency  of  cracks  as  impact  velocity  increases 
because  the  release  waves  initiated  by  failing  elements 
are  captured  by  the  underlying  finite  element  method. 
As  the  impact  velocity  is  increased,  the  release  waves, 
which  have  a  constant  velocity,  are  increasingly  incapa¬ 
ble  of  shielding  the  material  from  the  growth  of  addi¬ 
tional  cracks. 

The  Kayenta  model  in  ALEGRA  preserves  the 
completely  failed  material  in  the  mesh,  which  still  has 
strength,  and  can  be  recompressed.  Currently,  the  ma¬ 
terial  model  yield  surface  for  SiC  and  B4C  softens  down 
to  a  material  with  a  slope  similar  to  that  of  gravel  or 
sand.50  The  rate  of  the  limit  surface  collapse  is  based  on 
crack  speeds  and  the  corresponding  element  volumes 
such  that  the  energy  release  rate  is  relatively  mesh  in¬ 
sensitive. 

The  simulations  in  Fig.  1 1  represent  preliminary 
results  of  a  Kayenta  material  model  for  B4C.  The  JH 


model  for  B4C51  is  used  as  a  starting  point  for  defining 
the  median  limit  surface  to  be  used  in  the  Kayenta 
model.  The  Kayenta  model  requires  an  input  of  the 
time  it  takes  for  a  sample  of  the  reference  volume  to 
soften  completely.  It  is  assumed  that  these  times  are 
proportional  to  the  crack  growth  speed.  Approximate 
crack  speeds  determined  by  Ramesh  and  colleagues52,53 
are  used  for  both  SiC  and  B4C.  A  value  of  10%  of  the 
longitudinal  wave  speed  is  used  for  SiC,  and  3.5%  for 
B4C.  The  relationship  between  these  values  and  the  ac¬ 
tual  speeds  generated  by  the  model  is  yet  to  be  deter¬ 
mined,  and  future  revisions  of  the  model  may  allow 
direct  specification  of  the  crack-growth  speed. 

Sphere  indentation  experiments  with  B4C  exhibit  a 
greater  frequency  of  both  radial  cracking,  as  well  as  lat¬ 
eral  and  cone  cracking,  when  compared  with  SiC.  Using 
a  higher  Weibull  modulus  of  10,  along  with  the  lower 
crack  speed  mentioned  above  in  the  B4C  model  results 
in  a  good  match  to  the  experimental  results.  The  effect 
of  changing  the  Weibull  modulus  is  shown  in  Fig.  11. 
As  the  Weibull  modulus  increases,  the  variability  is  re¬ 
duced,  and  the  material  is  more  likely  to  fail  at  multiple 
locations  at  once  when  the  material  is  quickly  loaded  in 
the  impact  event. 


Fig.  11.  Weibull  modulus  effect  on  radial  cracking  in  boron  carbide  simulations  impacted  at  400  m/s. 
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This  indicates  that  it  may  be  possible  to  use  these 
types  of  experiments  to  calibrate  both  the  rate  param¬ 
eters  and  the  variability  parameters.  The  authors  hope  to 
obtain  characterization  experiments  of  B4C  similar  to 
those  completed  for  SiC  to  confirm  this  result. 

Conclusions 

In  this  paper,  dynamic  indentation  experiments 
have  been  assessed  for  their  usefulness  in  calibration 
and  validation  of  ceramic  material  models  that  account 
for  uncertainty  and  scale  effects  in  material  strength. 
Sphere  impact  experiments  between  50  and  500  m/s  are 
conducted,  which  promote  damage  and  failure  initia¬ 
tion,  without  substantial  penetration  or  complete  target 
failure.  Silicon  carbide  and  B4C  ceramic  targets  are 
analyzed  through  cross  sections  and  X-ray  computed 
tomography. 

The  number  of  radial  cracks  observed  in  the  dy¬ 
namic  indentation  offers  a  means  of  calibrating 
variability  in  strength  and  the  rate  of  the  progression 
of  damage.  Changes  in  the  character  of  cone  and  lateral 
cracking  that  occur  as  impact  velocity  increases  may 
serve  as  validation  for  damage  models,  provided  the 
experiments  are  not  used  for  the  initial  calibration. 

The  Kayenta  model  has  been  used  as  a  means  of 
exploring  nontraditional  features  in  modeling  dynamic 
failure  of  ceramics.  For  example,  third  invariant  depen¬ 
dence  allows  the  failure  criterion  to  transition  smoothly 
from  a  principal  stress  criterion  at  low  pressures  to 
capture  brittle  failure  behavior,  to  a  more  traditional 
von  Mises  plasticity  profile  at  high  pressures.  This  limit 
surface  profile  is  statistically  perturbed  to  reproduce 
observed  experimental  variability  in  strength  that  is 
pronounced  at  low  pressure  and  is  strongly  scale  depen¬ 
dent.  Softening  and  damage  is  modeled  by  allowing  the 
Kayenta  limit  surface  to  collapse  to  a  failed  surface  that 
is  representative  of  sand  or  gravel.  Calibrating  Kayenta’ s 
scale-sensitive  time-to-failure  criterion  requires  time- 
resolved  failure  data  to  supplement  the  dynamic  inden¬ 
tation  data. 

Mesh  sensitivity  continues  to  undermine  the  pre¬ 
dictive  capability  of  damage  models.  Incorporation  of 
scale  effects  and  uncertainty  dramatically  reduces  mesh 
sensitivity  in  dynamic  indentation.  These  features  also 
capture  trends  in  Brazilian  data,  but  do  not  eliminate 
mesh  dependence.  Accordingly,  advances  in  modeling 
and  validation  are  required  to  develop  material  models 


that  are  predictive  for  a  variety  of  loading  scenarios  us¬ 
ing  a  single  parameter  set. 
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